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The  materials  contained  in  this  report  were  presented  to  Dr.  Richard  G.  Brandt,  the 
sponsoring  Scientific  Officer,  as  part  of  an  October  9,  1991  site  visit  to  our  laboratory  in  the 
Physics  Department  at  Texas  A&M  University.  The  information  is  largely  pictorial  in  form, 
and  is  divided  into  four  sections.  Comments  pertinent  to  the  individual  figures  may  be  found 
in  each  of  Sections  II-IV. 

Section  I  begins  with  a  reprint  entitled  The  Image  Potential  in  Scanning  Tunneling 
Microscopy  of  Semiconductor  Surfaces  which  was  published  in  the  Journal  of  Vacuum 
Science  and  Technology  as  part  of  the  proceedings  for  PCSI-18.  A  comprehensive  and  more 
extensive  account  of  this  work  is  presented  in  an  accompanying  article  entitled  Internal  Image 
Potential  in  Semiconductors:  Effect  on  Scanning  Tunneling  Microscopy  which  was  recently 
submitted  to  Physical  Review.  These  papers  are  the  product  of  our  collaboration  with  the 
theoretical  group  of  Dr.  Roland  Allen. 

Section  II  traces  the  development  of  our  efforts  to  design  and  assemble  a  versatile  ultra 
high  vacuum  apparatus  for  multi-technique  surface  analysis,  tunneling  microscopy,  and  in-situ 
sample  preparation.  Aspects  of  the  design  process  facilitated  by  access  to  advanced  computer- 
based  software  packages  are  highlighted.  Particular  examples  include  three-dimensional 
visualization  of  the  complete  in-vacuo  sample  transfer  and  manipulation  scheme,  as  well  as 
finite-element  modeling  of  structural  modes  relevrnt  to  the  optimization  of  system  rigidity  for 
vibration  isolation. 

Test  data  demonstrating  operation  of  the  Perkin-Elmer  15-155  Cylindrical  Mirror 
Analyzer  installed  in  our  UHV  system  for  Auger  Electron  Spectroscopy  are  presented  in 
Section  III. 

Section  IV  documents  some  preliminary  imaging  and  spectroscopic  data  obtained  using  a 
Scanning  Tunneling  Microscope  (STM)  delivered  to  us  in  late  April  by  Omicron 
Vakuumphysik  .  This  instrument  has  been  specially  modified  to  facilitate  implementation  of 
Ballistic  Electron  Emission  Microscopy  (BEEM).  Photographs  detailing  the  flange-mounted 
STM  stage,  both  alone  and  in  relation  to  the  remainder  of  the  UHV  system,  may  be  found  in 
Section  I. 

The  individuals  responsible  for  carrying  out  our  ONR  sponsored  program  during  this 
past  year,  and  whose  efforts  have  led  to  the  results  reported  here,  are  listed  on  the  following 
page. 
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Section  I 


REPRINTS 


The  image  potential  in  scanning  tunneling  microscopy  of  semiconductor 
surfaces 

Z.-H.  Huang,  M.  Welmer.  and  B.  E.  Allen 

Department  of  Physics,  Texas  A^\f  University.  College  Station,  Texas  77843 
( Received  January  29,  1991;  accepted  April  17, 1991) 

A  semiclassical  screening  theory  for  semiconductors  is  used  to  treat  the  image  potential 
throughout  a  one-dimensional  metal-vacuum-semiconductor  junction  in  the  absence  of  surface 
states.  The  image  interaction  due  to  induced  surface  charge  at  the  semiconductor-vacuum 
interface  produces  an  effective  band  bending  in  the  semiconductor  interior,  whose  influence  on 
electron  tunneling  through  the  junction  is  investigated.  In  particular,  we  examine  the  effect  of  this 
extra  potential  on  STM  measurements  of  the  apparent  tunneling  barrier  at  unpinned 
semiconductor  surfaces,  and  consider  n-type,  H-terminated  Si(  1 1 1 )  as  a  specific  example.  In  the 
instance  where  tip-induced  band  bending  may  be  neglected,  we  find  the  role  of  the  image  potential 
in  metal-vacuum  semiconductor  tunneling  to  be  qualitatively  similar  to  its  role  in  metal-vacuum- 
metal  tunneling:  The  experimentally  determined  barrier  height— obtained  from  the  logarithmic 
derivative  of  the  tunneling  current  with  respect  to  tip-sample  separation — deviates  little  from  the 
electron  affinity  of  the  semiconductor,  while  the  theoretical  barrier  height  exhibits  a  stronger 
dependence  on  the  vacuum  gap  width.  The  origin  of  this  behavior  appears  to  lie  in  the  fact  that 
while  the  image  interaction  in  the  semiconductor  may  be  long  range,  the  image  induced  tunneling 
barrier  is  not. 


An  extra  electron  placed  in  a  polarizable  medium  always 
induces  a  redistribution  of  charge.  Near  a  surface  or  inter¬ 
face,  the  collective  response  of  other  electrons  to  this  exter¬ 
nal  perturbation  is  a  net  accumulation  of  charge  at  the  inter¬ 
face.  The  image  potential  describes  the  interaction  of  the 
extra  electron  with  that  induced  surface  charge. 

It  has  been  a  matter  of  some  interest  whether  or  not  the 
scanning  tunneling  microscope  (STM),  with  its  precisely 
controllable  vacuum  gap,  might  provide  a  unique  experi¬ 
mental  opportunity  to  study  the  role  of  image  interactions  in 
electron  tunneling.  The  question  has  been  thoroughly  ad¬ 
dressed  both  theoretically  and  experimentally  for  the  case  of 
mctal-vacuum-metal  (MVM)  tunneling.  In  the  simplest 
approximation,  one  models  the  STM  as  two  planar,  parallel, 
electrodes,  separated  by  a  distance  5,  and  calculates  the  clas¬ 
sical  correction  to  the  tunneling  potential  using  the  multiple- 
image  method.  *  ^  The  resulting  barrier  height  is  only  slightly 
affected  at  large  separations,  but  is  significantly  reduced 
when  5  is  the  order  of  a  few  angstroms.  Similar  conclusions 
may  be  drawn  from  more  sophisticated  density-functional 
calculations.  The  quantity  determined  in  STM  measure¬ 
ments.  however,  is  an  apparent  barrier  height,  <I>^ ,  defined 
by^ 

<*>^  =  ((  -  \/A)d\nI/ds]\  (1) 

Here  I  is  the  tunnel  current,  s  is  the  vacuum  gap  width,  and 
A  ^  1 .025  eV  “  A  *  ‘ .  Experimental  rcsults^’^  indicate  the 

logarithm  of  /  vs  5  is  nearly  linear  over  a  wide  range,  so  that 
is  essentially  independent  of  gap  width,  and  has  a  value 
close  to  the  average  work  function  of  tip  and  sample  elec¬ 
trodes.  As  Binnig  et  al.^  and  Coombs  et  ai^  explain,  it  ap¬ 
pears  that  <t>^  is  insensitive  to  the  image  interaction  because 
of  a  cancellation,  between  the  distance  dependent  barrier 
height  and  the  rate  of  change  of  this  barrier  height  with  in¬ 
creasing  gap  separation,  that  is  specific  to  the  form  of  the 
image  potential 


If  one  replaces  the  metal  sample  electrode  with  a  nonde¬ 
generate  semiconductor,  the  classical  image  force  in  the 
vacuum  is  slightly  reduced  due  to  the  semiconductor’s  finite 
polarizability,  but  this  does  not  qualitatively  alter  the  form 
of  the  potential  in  this  region.  A  more  important  distinction, 
however,  is  that  a  long  range  potential  also  extends  into  the 
semiconductor  interior,  and  this  results  in  an  effective  band 
bending  which  becomes  significant  for  electrons  near  the 
semiconductor-vacuum  interface.**’  This  band  bending  con¬ 
tributes  to  the  overall  tunneling  barrier  experienced  by  elec¬ 
trons  originating  from  the  bottom  of  the  conduction  band  m 
n-type  material.  It  has  previously  been  noted  that  the 
space  charge  induced  by  the  presence  of  a  metal  tip  near  an 
unpinned  semiconductor  surface  is  highly  sensitive  to  the 
distance  between  tip  and  sample,  and  that  varying  this  dis¬ 
tance  has  a  profound  effect  on  both  the  tunneling  probability 
and  the  experimentally  inferred  barrier  height.  We  may  in¬ 
quire,  then,  whether  the  effective  band  bending  associated 
with  image  interactions  in  a  metal-vacuum-semiconductor 
( MVS )  junction  is  similarly  sensitive  to  the  gap  spacing,  and 
therefore  presents  a  distinctive  signature  of  the  image  inter¬ 
action  which  is  unavailable  ir  MVM  tunneling.  It  is  the  pur¬ 
pose  of  this  paper  to  addrr  s  that  question. 

In  what  follows,  we  derive  a  semiclassical  expression  for 
the  image  potential  experienced  by  a  single  electron  travers¬ 
ing  an  ideal,  nondegcncrate  MVS  junction  without  surface 
states.  This  potential  is  used  to  calculate  a  tunneling  trans¬ 
mission  probability  by  solving  a  one  electron  Schrodinger 
equation  in  ,he  effective  mass  approximation.  An  apparent 
barrier  height  <I>^  is  then  extracted  with  the  aid  of  Eq.  (I )  for 
the  specific  case  of  an  n-type  hydrogen-terminated  Si(  1 11  )- 
Au  ystem  when  the  applied  bias  has  been  chosen  to  elimi¬ 
nate  tip-induced  band  bending. 

Expressions  for  the  classical  image  potential  in  an  ideal 
metal-insulator-semiconductor  (MIS)  junction  similar  to 
that  shown  in  Fig.  I  have  been  previously  denved It  is 
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semiconductor  vacuum  metal 


Fig  I  Energy  diagram  for  an  ideal  semiconductor-vacuum-meul  junction, 
when  the  applied  bias  V  has  been  chosen  to  eliminate  any  electric  field 
ansing  from  the  difference  in  metal  and  semiconductor  work  functions, 
^  d, )  Here  x  the  semiconductor  electron  affinity,  and  d,  the 
offset  of  the  Fermi  level  from  the  bottom  of  the  conduction  band. 

weJI-known,  however,  that  this  potential  diverges  at  the  in¬ 
terfaces.  The  unphysical  divergences  may  be  circumvented 
using  a  Thomas-Fermi  theory  of  dielectric  screening.' ‘ 
The  Thomas-Fermi  dielectric  function  in  a  bulk  metal  is 
e^ik)  =  \  -h  yi,/k  where  is  the  corresponding  Thom¬ 
as-Fermi  wave  vector.  For  a  point  charge  located  at  posi¬ 
tion  To ,  the  electrostatic  potential  inside  the  metal,  { r.Fo ) , 

satisfies  the  equation 

=  -4vqg6{T -To).  (2) 

Within  the  same  approximation,  the  screening  response  of  a 
bulk  semiconductor  to  a  point  charge  may  be  modeled  using 
the  dielectric  function  given  by  Inkson," 

e,(k)  =  1  4-  [(f,  -  1)/{1  4-6,/c  Vy^)],  (3) 

where  is  the  static  dielectric  constant  and  y,  an  effective 
Thomas-Fermi  wave  vector.  It  is  easily  shown  that  the  elec¬ 
trostatic  potential  inside  the  semiconductor,  ^  r,ro ),  then 
satisfies  the  equation ' ' 

(V^  -  }i)ii,(T,ro) 

=  -  (}i/e,)u(T,To)  -  ^tr(qo/€,)S(r  -To),  (4) 


inside  the  semic^^nductor,  and  then  determine  6^  from  Eq. 
(4).  With  this  choice  for  u,  the  resulting  semiclassical  poten¬ 
tial,  d,,  naturally  approaches  the  classical  solution  at  dis¬ 
tances  greater  than  a  few  Thomas-Fermi  wavelengths  from 
the  interface.  To  obtain  the  semiclassical  image  potential, 
one  simply  evaluates  the  semiclassical  electrostatic  potential 
at  the  position  of  the  point  charge,  excluding  the  field  pro¬ 
duced  by  the  charge  itself. 

Defining  the  quantities 


a,(k) 


k 

k 

ik-^fA'  '-  ' 

1  -an,(k) 

1  +a„(k)  ' 
1  -aAk) 

1  T  a,{k) 


I  5a) 


i5b) 


( 6a ) 


I  6b) 


we  show  in  the  Appendix  that  the  semiclassical  image  poten¬ 
tial  felt  by  a  point  charge  in  the  vacuum  region,  0  <  r,,  <  s, 
is 


<Pc 

where 


;,(Zo)=9of  (Bte'^'  +  C^e^''“')dk. 

Jo 


,  ri,F  ./€,  —  77. e  e 

Bl  =  - 

e- 

r’v  /  ic  (  f  -  4,  J  ,  -  4  (  j  -  i ,  . 

^  F\/e,  -  e  V„e _ 


Cl  =  T]^e‘ 
and 

with 


-  Z.A 

e  —  e 


\  -0e 


(1) 


(8) 

{^) 


(  10) 


(£,-!)/(£,  + 1). 


(ID 


where  u(r,ro )  is  the  vacuum  potential,  go/\T  —  Tq\.  Equa¬ 
tions  (2)  and  (4)  are  just  coordinate  space  representations 
of  the  fundamental  linear  relationship, 

=  u{k)/€^Jk). 

In  turning  to  the  junction  shown  in  Fig.  1,  we  assume  the 
dielectnc  functions  defined  for  bulk  media  are  unaffected  by 
the  presence  of  any  interfaces,  so  that  Eqs.  (2)  and  (4)  re¬ 
main  valid  throughout  their  respective  medium.  The  semi¬ 
classical  solution  to  the  electrostatic  potential  inside  the  met¬ 
al.  (^>^,  is  then  given  by  Eq.  (2),  subject  to  boundary 
conditions  at  the  metal-vacuum  interface.  In  order  to  deter¬ 
mine  the  semiclassical  expression  for  however,  we  must 
specify  the  vacuum  potential,  u,  in  addition  to  boundary 
conditions  at  the  semiconductor-vacuum  interface.  In  evalu¬ 
ating  u,  we  need  to  consider  the  contribution  of  all  charge 
external  to  the  semiconductor  half  space,  in  addition  to  any 
charge  placed  inside  the  medium.  As  a  lowest  order  approxi¬ 
mation,  we  add  to  u  the  vacuum  potential  from  the  external 
image  charges  which  generate  the  correct  classical  solution 


The  corresponding  potential  energy,  (Zo  ),  is  then  giv¬ 
en  by 

The  image  potential  for  semiconductor  ( z,,  <  0 ) 

is  similarly  found  to  be 


<PA2o) 


—  f  {F\e'^  ^  ^  y:)^  ^^z^]}dk. 


where 

F\  = 
and 


-  -  p 

^  -  p 


(13) 


Al  =  (2a,/{l  4-  a,){r)^f  i }  +  (f,  -  1) 

^  -  V,  ]  exp[(<c-  +  rf)'  •;„]  )/ 

“3  (14) 
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Finally  for^,,  in  the  metal  ( j,,  >  o  we  have 

=  q.,  I  e.xp  [  (  A:  -  -  ) '  ' ( s  -  j,, )  ]  Z)  '”dk. 

(15) 

With 

D:=  ,-e.xp[(<l:^-;i)(5-z„)).  (16) 

1  -  7.7^^ 

It  IS  easily  venfied  that  the  expressions  above  reproduce 
the  classical  image  potential  in  the  limits  —  x.  If  we 
consider  instead  the  limit  of  very  large  gap  spacing,  5  —  x, 
we  recover  the  semiclassical  result  of  Newns'*  for  an  isolat¬ 
ed  vacuum-metal  interface,  and  an  analogous  result  for  the 
semiconductor-vacuum  interface.  In  the  limit  of  a  vanishing 
vacuum  gap.  5  —  0.  our  junction  reduces  to  a  semiconductor- 
metal  contact.  In  this  case,  the  results  derived  above  differ  in 
detail  from  those  quoted  by  Inkson'*  because  of  our  differ¬ 
ent  ansatz  for  w(  r.r^  ).  In  effect.  Inkson  treats  the  semicon¬ 
ductor-metal  interface  by  constructing  u  from  the  solution 
for  a  classical  dielectric  in  contact  with  a  Thomas-Fermi 
metal.  Our  starting  point  for  the  same  problem  would  be  the 
solution  to  a  classical  dielectric  in  contact  with  a  classical 
metal.  If  we  mtroduce  a  vacuum  gap  betw'een  the  two  media, 
the  distinction  between  these  approaches  becomes  unimpor¬ 
tant  for  separations  s  which  are  large  compared  with  a 
Thomas-Fermi  screening  length  in  the  metal. 

The  effective  one-electron  potential  in  the  junction  is  ob¬ 
tained  by  adding  the  semiclassical  image  potential  to  the 
band  diagram  shown  in  Fig.  1 .  To  clearly  separate  the  image 
interaction  from  other  factors  w  hich  might  influence  the  de¬ 
termination  of  an  apparent  barner  height, we  presume  the 
applied  bias  V  has  been  chosen  to  eliminate  the  electric  field 
between  tip  and  sample  in  the  absence  of  image  corrections. 
Selecting  the  bottom  of  the  semiconductor  conduction  band 
as  our  energy  reference,  the  one-electron  potential  is  then 
given  by  \qo<f,  m  the  semiconductor,  X  vacu¬ 

um,  ana  \qn(frr,  ^  -h  £'^)  in  the  metal.  This  func¬ 
tion,  though  finite,  is  discontinuous  at  the  semiconductor- 
vacuum  and  vacuum-metal  boundaries.  It  is  unreasonable  to 
expect  our  semiclassical  treatment  to  be  entirely  accurate 
within  a  screening  length  of  the  interfaces  for  a  variety  of 
reasons.  First,  the  treatment  assumes  bulk  values  for  the 
screening  parameters.  Second,  the  correct  behavior  of  the 
image  potential  across  an  interface  can  only  be  deduced  by 
taking  full  account  of  many-body  effects.*^  Finally,  there  is 
no  a  prion  relationship  between  the  effective  band  bending 
we  calculate  in  the  immediate  vicinity  of  the  interfaces,  and 
the  one-electron  potential  we  use  to  model  band  structure, 
which  guarantees  that  the  contribution  from  the  surface  di¬ 
pole  layer,  together  with  exchange  and  correlation  effects, 
will  produce  a  smoothly  varying  total  potential  across  the 
interface.  Since  we  believe  the  effective  one-electron  poten¬ 
tial  should  be  continuous  when  all  these  effects  are  self-con- 
sistently  accounted  for.  we  have  supplemented  our  picture 
by  employing  a  third-order  polynomial  to  smoothly  join  the 
potential  and  its  denvative  at  points  a  distance  from 
either  side  of  the  boundaries.  We  then  use  this  interpolation 
to  represent  the  tunneling  potential  within  the  interfacial 


Fic.  2.  Effective  one-electron  potential  energy  as  a  function  of  position 
Solid  line:  including  the  semiclassical  image  potential  with  interpolation 
Dotted  line:  including  the  classical  image  interaction.  Dashed  line  neglect¬ 
ing  image  effects.  Here  ^  =  43  cV,  =  5. 1  cV,  and  £^  =  53  eV 

layer.  An  example  of  this  interpolated  potential  is  presented 
in  Fig.  2.  where  we  contrast  our  results  with  the  classical 
calculation,  for  a  junction  consisting  of  a  ^-type  hydrogen- 
terminated  Si(lll)  surface  and  gold  tip  electrode,  wuh 
y^  =  1.7  A- =  1.16  A  ',  and  F,  -  11.9.'**^'^ 

The  method  used  in  computing  the  tunneling  current  is 
the  same  as  that  described  in  Ref  10.  Bnefly.  the  potential  is 
substituted  into  a  one-dimensional,  effective  mass  Schro- 
dinger  equation,  which  is  solved  numerically  to  evaluate  a 
transmission  probability.  Finite  temperature  supply  func¬ 
tions  are  then  employed  to  determine  the  current  density  ’ 
The  apparent  barrier  as  a  function  of  gap  spacing  is  calculat¬ 
ed  via  Eq.  (1).  and  our  results  illustrated  in  Fig.  3  The 
theoretical  barrier  height,  defined  by  the  maximum  of  the 
potential  in  the  vacuum  region,  is  also  shown  for  compan- 
son.  As  one  sees  from  the  figure,  the  image  interaction  re¬ 
duces  the  theoretical  barrier  height  with  decreasing  gap 
spacing,  but  the  effect  is  significant  only  for  5  less  than  10  A. 
The  apparent  barrier  height  defined  in  Eq.  (I ),  on  the  other 


Fig  3  Apparent  {dotted  line)  and  theoretical  (solid  line)  barner  height  vs 
electrode  separation  r  for  T  -  300  K  The  results  are  vinually  independent 
of  doping  in  the  nondegenerate  regime 
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Fig.  4  Differential  current  density  disinbution  vs  normal  energy  at  several 
tip-sample  separations 

hand,  never  vanes  appreciably,  and  differs  little  from  the 
electron  affinity  of  the  semiconductor.  This  is  essentially 
equivalent  to  the  result  for  a  MVM  tunneling  junction. 

The  similanty  in  apparent  barrier  height  as  a  function  of 
gap  spacing  in  MVS  and  MVM  junctions  suggests  that  the 
influence  of  image  interactions  on  electron  tunneling  in  the 
semiconductor  is  not  significant.  This  appears  contrary  to 
expectations  based  on  the  long-range  nature  of  the  image 
force  m  a  dielectric  medium.  Specifically,  since  tunneling 
electrons  in  a  nondegenerate,  n-type  semiconductor  origi¬ 
nate  from  the  bottom  of  the  conduction  band  ( see,  for  exam¬ 
ple,  the  calculated  differential  current  density  distributions 
shown  in  Fig.  4 ) ,  one  would  argue  that  the  tunneling  current 
should  be  sensitive  to  vanations  in  the  precise  form  of  this 
long-range  potential  with  electrode  separation.  However, 
there  is  an  interesting  point  whose  importance  is  overlooked 
in  prior  classical  treatments  of  this  problem,®-^  but  which  is 
crucial  to  the  discussion  here:  Near  a  semiconductor  surface, 
the  image  potential  is  dominated  by  induced  charge  at  the 
semiconductor-vacuum  interface.  The  sign  of  this  charge  de¬ 
pends  only  on  the  discontinuity  in  polarizability  at  the  inter¬ 
face  and  results  in  a  repulsive  force.  At  distances  in  the  semi¬ 
conductor  intenor  which  arc  large  compared  to  the 
electrode  spacing,  however,  the  vacuum-metal  interface  also 
contributes  to  the  induced  charge  determining  the  long- 
range  potential.  Since  the  polarizability  of  the  metal  exceeds 
that  of  the  semiconductor,  the  net  force  is  now  attractive. 
Because  the  image  force  must  vanish  far  from  the  semicon¬ 
ductor  surface,  at  distances  between  these  two  extremes  the 
image  potential,  and  therefore  the  effective  band  bending, 
must  change  sign.  Consequently,  the  induced  tunneling  bar- 
ner  only  extends  to  the  point  where  the  image  potential 
passes  through  zero,  and  is  thus  of  finite  range.  The  effect  of 
the  remaining  long-range  attractive  part  of  the  image  poten¬ 
tial  on  electron  tunneling  may  be  seen  in  the  resonance  struc¬ 
ture  displayed  by  the  differential  current  densities  shown  in 
Fig.  4  In  the  particular  case  considered  here,  the  zero  cross¬ 
ing  IS  given  by  =  3.59.  For  all  of  our  calculations,  we 
have  truncated  the  effective  one-dimensional  potential  on 
the  semiconductor  side  of  the  junction  at  a  distance  which  is 
large  compared  with  Zo  •  ^i^d  presumed  that  Debye-Huckel 


screening  due  to  free  earners  is  insignificant  on  this  length 
scale. 

In  conclusion,  we  have  examined  a  model  STM  junction 
in  which  the  only  contnbution  to  the  tunneling  barner  in  the 
semiconductor  is  from  the  image  potential.  In  this  case  the 
apparent  barrier  height  is  insensitive  to  electrode  separation. 
The  finite  range  of  the  image  induced  barner  in  the  semicon¬ 
ductor  provides  a  natural  explanation  for  this.  As  seen  m 
Fig.  2,  the  total  area  contnbuted  by  the  tunneling  barner  in 
the  semiconductor  is  typically  small  compared  wuh  that 
contributed  by  the  vacuum  gap.  Furthermore,  the  small  ef¬ 
fective  mass  for  electron  tunneling  along  the  (111)  axes  in 
silicon  enhances  the  transparency  of  any  barner  in  this  direc¬ 
tion. 

Despite  fundamental  differences  in  the  long-range  screen¬ 
ing  propenies  of  metals  and  semiconductors  which  suggest  a 
more  prominent  role  for  the  image  potential  in  metal-vacu- 
um-semiconductor  versus  metal-vacuum-metal  tunneling,  it 
appears  that  the  experimentally  determined  barner  height  in 
STM  measurements  will  be  equally  insensitive  to  image  ef¬ 
fects  in  both  instances.  The  physical  origin  of  this  insensiti¬ 
vity  is  different  in  the  two  cases,  however,  and,  for  the  semi¬ 
conductor,  may  be  traced  to  the  short-range  nature  of  the 
tunneling  barrier  induced  by  the  long-range  image  force  in  a 
dielectric  medium  near  a  metal  electrode. 

We  thank  J.  Kramar  for  providing  pan  of  the  computer 
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of  Naval  Research  through  Grant  Nos.  N00014-89-J-3029 
and  N00014-91-J-1 126.  and  by  the  Robert  A.  Welch  Foun¬ 
dation. 

APPENDIX 

We  first  consider  the  case  where  a  point  charge  is 
placed  in  the  vacuum  region  of  the  junction  shown  in  Fig.  1 
The  electrostatic  potential  in  each  of  the  three  regions  is 
determined  by  the  equations 

-f,)6,{r,To)  =  -  (}j/€,)u{T.T,y), 

z<Q  { semiconductor ) 

(r.ro  )  =  -  ^vq^SiT  -r^), 

*  ( .A  I  ) 

0<z<s  (vacuum ) 

(V^  -  yi)<#„(r,ro)  =0, 

^  z>5  (metal) 

where  Tq  =  (0,0,Zo )  is  the  position  of  q^.  The  general  solu¬ 
tion  to  these  equations  subject  to  the  boundary  condition 
that  both  and  6^  vanish  at  x .  may  be  expressed  in  cylin¬ 
drical  coordinates  as  follows; 

€,  Jo  6, 

xf  i/)(/cp)  exp[  (^  ^ -h  K)‘  ‘2]dA.  (  .A2) 
Jo 

(r.ro)=^o  f 

Jn 

^C\e  ^‘)dk.  iA3) 
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and 

r  * 

(r,ro)=9o  D[J„{>(p)  exp[  -  {k-  -r  yi,)'  'zjdk. 

JO 

(A4) 

The  first  term  in  Eq.  (A2)  is  the  particular  solution 
wlr.fo  )/€,,  where  we  have  approximated  u  using  the  well- 
known  classical  expression"*  given  in  Eq.  (10).  The  second 
term  is  a  homogeneous  solution  ansing  from  short  range 
screening  in  the  semiconductor.  Equation  (A3)  forc^,,  con¬ 
tains  a  particular  solution  due  to  the  point  charge  qQ,  as  well 
as  a  homogeneous  solution.  The  coefficients  A  I,  Bl, 
and  D I  are  determined  by  four  boundary  conditions  at  the 
interfaces,  namely 


d(i>^ 

dz 


dz 


(A5) 

(A6) 

(A7) 

(A8) 


d6,  j 

17 1. - 

and 

36  j 

In  imposing  continuity  of  the  electrostatic  potential  and  its 
denvative  at  these  boundanes,  we  make  the  physical  as¬ 
sumption  that  all  interfacial  charge  is  distributed  over  a  fi¬ 
nite  width.  Since  we  exclude  an  abrupt  surface  charge,  the 
factor  €  which  normally  accompanies  the  classical  continu¬ 
ity  equation  for  d<b/dz  no  longer  appears  in  our  boundary 
conditions  for  the  semiclassical  solution. 

The  image  potential  is  obtained  from  6^.  by  setting  r  = 
and  excluding  the  self-energy  term.  The  resulting  expression 
IS  given  by  Eqs.  ( 7)--(  1 1 ). 

When  ^0  is  in  the  semiconductor,  the  regional  differential 
equations  become 

f  (V-  -  r.r„  )  =  -  (K/e,  )tt(r.ro) 

-  4n-<f^S(r  -  To). 
z<0  (semiconductor); 

(r,ro  )  =  0, 

0<7<s  (vacuum); 

(V’  -  }i,)(i^(r,ra)  =0, 

.  z>s  ( metal ) 


(A9) 


The  solutions  in  these  three  regions  are 
6  i  rr,) )  -  ^  r  Jo(^p)e^  '  -h  qo 

€,  Jo 


€.  -  I 


[  Jn(kp)  exp[(^‘  -  z,^:]dk 

Jo 

-f 

€,  Jo 


F[J„(lcp)e‘-^dk  +  — 


^  * 

X  A  lJ:){kp)  exp[(/c'  J-  ^z]dk, 
Jo 


(AlO) 
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Jo 

(A12) 

Additional  inhomogeneous  terms  now  arise  in  Eq.  (AlO) 
from  the  presence  of  the  point  charge  in  the  semiconduc¬ 
tor.  “  The  classical  expression"*  *’  is  again  used  to  represent 
the  vacuum  potential  u  and  the  resulting  image  potential  is 
given  by  Eqs.  ( I2)-(  14). 

Finally,  when  q^  is  m  the  metal,  we  have 

''(V-  -  )  =  -  (r':/e,  }w(r,ro  ), 

r<0  (semiconductor); 

0<z<5  (vacuum); 
iV"  -  ri)d>^{r,ro)  =  -  47rq,,6ir  - 
^  z>s  { metal ) . 

whose  general  solutions  are 

dArSo)  =  —  f*  F:J„{kp)e^^dk  ^  — 

6,  Jo  e 

^  f  ^  >7) '  '^]3k. 

Jo 

(A14) 

<^..(r,ro) -^0  f  (A15) 

Jo 

and 

d>m0.rn)=go  I  Joikp) 

Jo 

X  exp [  —  ( '  2  -  Zo  ]dk 
^0  f  ^  T^o  ( kp) 

Jo 

xexp[  -  ik^^  -t  yi,)^  'z]dk.  (Alb) 

Here  we  have  assumed  the  vacuum  potential  u  vanishes  in 
the  semiconductor,  as  it  must  in  the  classical  case  where  the 
metal  completely  screens  the  point  charge.  Similar  deriva¬ 
tions  then  result  in  Eqs.  (15)  and  (16) 


’  J  G.  Simmons.  J  Appl  Phys.  34,  1793  {  1963) 

Q.  Nguyen,  T-  E.  Fcuchiwang,  and  P  H.  Cutler.  J  Phys  47.  C2 
(1986). 

'N  D  Lang.  Phys.  Rev  B37.  10395  ( 1988) 

*G.  Binmg,  N  Garcia.  H,  Rohrer.  J.  M.  Soler.  and  F  Flores.  Phys  Rev  B 
30.  4816  (1984) 

’G.  Binnig,  H.  Rohrer.  Ch.  Gerber,  and  E.  Weibcl,  Appl.  Phys.  Lett  40. 
178  (1982). 

M  H.  Coombs.  M  E.  Welland,  and  J  B  Pcthica,  Surf  Sci  198.  L353 
(1988) 

^  For  treatments  in  nonplanar  geometnes  see.  for  example,  A  A  Lucas 
ef  al,  }  Phys.  45.  C9  ( 1984);  H.  Morawiiz  et  a/.  Surf.  Sci  180.  333 
(  1987);  J  Mahanty  and  M  T  MichaJewter,  Austral  J  Phys  40,  413 
(1987). 

^  M.  Kleefstra  and  G  C  Herman.  J  Appl.  Phys.  51.  4923  ( 1980) 

A  Tuguiea  and  D  Dascalu.  J  Appl.  Phys  56.  2823  (  1984) 

M.  Wcimer.  J  Kramar.  and  J  D  Baldcschwicler,  Phys  Rev  B  39.  55^2 
( 1989). 


2404  Huang,  Weimer,  and  Allen:  STM  of  semiconductor  surfaces 

J  C  InkM'n.  J  Ph>s  C4. 

D  Vt  J  Chcrn  Ph\s  50.  45'2  ' 

D  Lang,  Solid  Slate  Hhvs  28.  , 

*  *S  W  A>hcrort  and  \  D  Mermin.  SVd.v  v/to  i  Holt.  Nesv  ^’ork. 
Id'M 

‘J  P  Waiter  and  M  L  Cohen.  Phys  Rev  B  2.  is:i  i  Id^O) 


2404 

"S  M  Sze.  Physics  of  SemiK.onducior  Deuces  \  W  iley.  New  '*  .■>rk.  !'■>/' 1 
The  doping  density  n  e'.seniialiy  irreiesanl  to  the  caicuiatim  presented 
here,  since  the  only  quantity  which  depends  on  n  is  the  povi’:.-n  i  :nc 
semiconductor  Fermi  level  In  the  nonde^eneraie  case  this  mereiv  atfevt-* 
the  overall  normalization  of  the  supply  functions 


J.  Vac.  Sci.  TechnoL  B.  Vol.  9,  No.  4,  Jul/Aug  1991 


Internal  Image  Potential  in  Semiconductors: 
Effect  on  Scanning  Tunneling  Microscopy 
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Abstract 

The  tunneling  of  electrons  from  a  semiconductor  surface  to  a  metal  tip,  across  a  vacuum  gap,  is 
influenced  by  two  image  interactions:  an  attractive  image  potential  in  the  vacuum  region,  which 
lowers  the  apparent  tunneling  barrier,  and  a  repulsive  image  potential  in  the  semiconductor  interior, 
which  raises  it  for  conduction  band  elections.  We  report  on  detailed  calculations  of  tunneling 
currents  and  apparent  barrier  heights  for  a  model  metal-vacuum-semiconductor  junction  which 
utilizes  semiclassical  dielectric  functions  to  compute  the  image  potential  in  all  three  regions.  The 
effect  of  image  forces  is  found  to  be  small  compared  to  that  of  either  the  vacuum  barrier  or  tip- 
induced  band  bending.  In  particular,  the  image-induced  barrier  in  the  semiconductor  has  only  a 
minor  itifluence  on  either  the  apparent  barrier  height  or  the  shape  of  current- voltage  characteristics, 
both  of  which  are  routinely  measured  in  STM  experiments.  This  finding  is  explained  by  a 
qualitative  WKB  analysis  and  several  simple  arguments. 
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I.  Introduction 


Electron  tunneling  through  a  junction  is  often  treated  by  considering  independent  electrons  moving 
in  an  effective  potential.  The  simplest  model  potential  for  a  metal-vacuum-metal  (MVM)  junction  is 
the  rectangular  barrier,  but  it  is  well  known  that  image  forces  modify  this  potential  significantly. 
Calculations  based  on  classical  electrostatics,  for  example,  reveal  that  the  image  interaction  reduces 
the  effective  vacuum  barrier  by  lowering  its  maximum  height  and  slightly  narrowing  its  width 
More  realistic  quantum  mechanical  calculations  employing  the  density  lunctional  formalism  indicate 
that  the  lowest  order  correction  to  this  classical  picture  is  simply  a  displacement  of  the  effective 
image  plane  by  several  tenths  of  an  angstrom,  along  with  a  smooth  matching  of  the  vacuum 
potential  outside  the  surface  to  the  minimum  of  the  conduction  band  in  the  me*''^  interior  The 
resulting  vacuum  barrier  height  now  depends  on  the  gap  spacing  s,  and  is  well  approximated  by 

=  (I) 

where  0o  represents  the  average  work  function  of  the  metal  electrodes,  So  accounts  for  the  shift  in 
image  plane  positions,  and  a~  10  eV  A. 

The  scanning  tunneling  microscope  (STM),  with  which  one  can  precisely  adjust  s,  and  hence 
continuously  vary  the  vacuum  barrier,  is  a  potentially  valuable  tool  for  investigating  this  effect  By 
modulating  the  distance  between  tip  and  sample  at  a  frequency  outside  the  microscope's  feedback 
bandwidth,  and  detecting  a  synchronous  ac  tunnel  current,  one  can  infer  an  apparent  barrier  height 
according  to 

<t>,  =[{-{/  A)d\n I  /dsf,  (2) 

where  -4  =1.025  eV'^/^A'l.  The  utility  of  (2)  stems  both  from  the  fact  that  it  is  an  experimentally 
accessible  quantity,  and  that  in  a  Wentzel-Kramers-Brillouin  (WKB)  approximation  for  the 
rectangular  barrier  it  yields  the  correct  barrier  height.  Early  experiments  on  metals,  however, 
revealed  no  discernible  reduction  in  by  image  forces  at  separations  of  more  than  a  few 
angstroms.  This  insensitivity,  it  was  argued  is  due  not  to  the  absence  of  image  effects,  but, 
rather,  to  the  particular  form  of  the  potential  in  (1),  for  which  a  WKB  analysis  shows  that  the  first 
order  term  in  an  expansion  of  <Pa  in  powers  of  5'^  vanishes.  Other  lines  of  reasoning  similarly 
supported  the  conclusion  that  the  apparent  barrier  height  (2)  is  unlikely  to  be  a  sensitive  probe  of 
the  image  potential  in  M'VM  tunneling  at  large  distances  Nevertheless,  the  question  of  image 
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effects  and  their  relevance  to  STM  experiments  continues  to  hold  great  interest.  Subsequent 
investigators  have  addressed  a  number  of  fundamental  issues,  including  the  role  of  non-planar 
geometry  the  precise  nature  of  the  effective  vacuum  barrier  at  short  distances 

the  influence  of  the  dynamic  image  interaction  in  time-dependent  tunneling 


If  we  replace  one  of  the  metal  electrodes  by  a  semiconductor  to  form  a  metal-vacuum- 
semiconductor  (MVS)  junction,  then  the  image  correction  to  the  vacuum  potential  at  large  distances 
will  be  similar  to  that  of  the  MVM  junction  when  scaled  to  reflect  the  dielectric  behavior  of  the 
semiconductor  A  qualitatively  new  phenomenon  now  appears,  however,  in  that  tunneling 
electrons  experience  an  additional  image  force  in  the  semiconductor  interior  due  to  the  finite 
polarizability  of  that  medium.  In  a  classical  multiple-image  analysis  of  the  planar  MVS  or  SVS 
junction,  tills  "intemal"  image  potential  is  given  by  the  expression 


Kiz) 


- ^ ^ 

[ns-z  (n  +  Ds-z) 


(3) 


for  an  extra  point  charge  q  located  at  z  <  0,  with  z  =  0  the  position  of  the  semiconductor-vacuum 
interface,  and  z  =s  the  position  of  the  vacuum-metal  interface.  Here,  £s  is  the  static  dielectric 
permittivity  of  the  semiconductor  and 


P  = 


i£s  +  £o)' 


(4) 


with  a  corresponding  expression  for  P'.  For  an  MVS  junction,  /?'=  1.  The  intemal  image 
potential  (3)  leads  to  an  effective  increase  in  the  semiconductor  band  gap  near  an  ideal 
semiconductor-vacuum  interface,  since  electrons  in  the  conduction  band  are  repelled  from  the 
surface  whereas  electrons  in  the  valence  band  are  attracted  Image-induced  features  of  this 
kind  have,  in  fact,  recently  been  proposed  as  an  explanation  for  the  thickness-dependent  energy 
shifts  observed  in  both  XPS  and  AES  spectra  of  thin  Si02  fdms  on  silicon  substrates  The 
situation  is  exactly  reversed  at  an  ideal  semiconductor-metal  interface  where  the  band  gap  is 
effectively  decreased,  a  phenomenon  already  anticipated  from  the  quantum  many-body  point  of 
view 


To  explicitly  show  the  functional  dependence  on  gap  spacing,  equation  (3)  may  be  written  in  scaled 
formas)  as 


2 


(5) 


where  q=zls  and  Us  is  a  dimensionless  function  independent  of  s.  In  an  MVS  junction,  V siU 
changes  sign  at  a  position  Zq  =  s^o(P)  because  the  contribution  from  polarization  charge  at  the 
metal-vacuum  interface  eventually  dominates  that  of  the  semiconductor-vacuum  interface  deep  in 
the  semiconductor  interior.  Thus,  a  conduction  band  electron  approaching  the  surface  will 
experience  an  image  interaction  which  is  at  first  attractive,  and  then  repulsive,  with 

(6) 

U,>0,^>l(P) 


The  image  induced  tunneling  barrier  in  the  semiconductor  is  then  strongly  dependent  on  electrode 
separation.  It  is  appealing,  therefore,  to  examine  the  possibility  of  detecting  this  additional  barrier 
using  the  scanning  tunneling  microscope. 

The  classical  expression  for  the  image  potential  is  inconvenient  to  use  in  tunneling  calculations 
implementing  transfer  matrix  methods,  since  it  diverges  at  the  interfaces,  where  unphysical  bound 
states  may  be  introduced.  A  semiclassical  theory  for  the  MVS  junction  which  removes  these 
divergences,  but  nevertheless  requires  an  interpolation  for  the  total  effective  potential  in  the  near 
interfacial  region,  has  recently  been  discussed  Initial  calculations  based  on  this  theory  showed 
that,  for  the  special  case  of  vanishing  electric  field  between  tip  and  sample,  the  contribution  to  0A 
from  the  image  induced  barrier  in  the  semiconductor  is  small.  The  insensitivity  of  0A  to  this 
additional  correction,  it  was  argued,  could  be  ascribed  to  the  short-range  nature  of  the  induced 
barrier  in  the  semiconductor,  which  therefore  appears  nearly  transparent.  It  is  the  purpose  of  this 
paper  to  re-examine  the  question  of  observing  such  image  effects  in  an  MVS  junction  under  more 
general  circumstances. 

To  do  so,  we  first  review  briefly  the  phenomenology  of  apparent  barrier  heights  in  tunneling 
experiments  at  unpinned  semiconductor  surfaces.  At  zero  applied  bias  there  will,  in  general,  be  a 
space  charge  induced  in  the  semiconductor  by  the  difference  between  tip  and  sample  work 
functions,  A<t>  =  as  shown  in  Fig  1(a)  for  n-type  material  in  depletion.  This  space 

charge  may  be  eliminated  only  by  applying  an  appropriate  bias  voltage,  Vpg  =  to  recover  the 
flat-band  condition  illustrated  in  Fig.  1(b).  Tip-induced  band  bending  of  the  kind  indicated  in  Fig. 
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1(a)  will  dramatically  affect  STM  measurements  of  0^,  for  reasons  which  are  easy  to  appreciate 
The  surface  potential  Vd  is  a  monotonically  decreasing  function  of  s,  with  a  maximum  at  5  = 
0.  Since  the  depletion  width  W  is  typically  much  longer  than  s,  only  those  carriers  with  a 
thermally  activated  kinetic  energy  greater  than  Vd  will  surmount  the  diffusion  barrier  to  arrive  at  the 
semiconductor-vacuum  interface.  Increasing  the  separation  between  tip  and  sample  exponentially 
suppresses  tunneling  through  the  vacuum  barrier  as  before,  but  now  also  exponentially  enhances 
the  number  of  electrons  reaching  the  semiconductor  surface.  Because  of  these  two  competing 
factors,  the  observed  variation  in  current  density  is  less  than  expected;  this  is  then  interpreted  as  a 
smaller  barrier  height  in  Eq.  (2).  If  the  image-induced  potential  (5)  is  superimposed  on  the  space 
charge  potential  of  Fig  1(a),  it  too  will  change  with  gap  spacing.  Since  the  energy  scale 
characterizing  the  carrier  distribution  is  kT,  we  may  now  ask  whether  inclusion  of  the 
semiconductor  image  potential  further  perturbs  the  carrier  distribution,  and  thus  the  logarithmic 
derivative  of  the  tunneling  current,  in  any  measurable  way. 

In  what  follows,  we  first  present  key  elements  of  the  semiclassical  image  potential  theory,  and 
describe  our  method  of  calculation,  in  Section  11.  Predictions  for  the  apparent  barrier  height  and 
tunneling  current-voltage  characteristics,  including  image  effects,  are  presented  in  Section  III  for 
the  particular  case  of  a  planar,  n-type  hydrogen-passivated  Si  (lll)/Au  junction.  We  devote 
Section  IV  to  a  qualitative  discussion  of  these  results  within  the  framework  of  a  WKB  analysis, 
which  permits  us  to  generalize  our  findings  beyond  the  specific  materials  we  have  chosen  to  study 
in  detail.  Section  V  summarizes  our  conclusions. 

II.  Description  of  the  Calculation 

Our  calculations  are  carried  out  under  the  assumption  that  the  system  consists  of  independent 
electrons  moving  in  an  effective  one-dimensional  potential.  The  bulk  metal  is  treated  as  a  free- 
electron  Fermi  gas  characterized  by  Fermi  energy  Ep  and  work  function  0.  The  semiconductor  is 
modeled  in  a  one-band  effective-mass  approximation  including  anisotropy.  When  assembled  into 
an  MVS  Junction,  the  semiconductor-vacuum  and  vacuum-metal  interfaces  perturb  these  bulk 
effective  potentials;  the  corrections  we  specifically  consider  here  are  introduced  by  the  multiple- 
image  interaction  and  tip-induced  space  charges.  Modifications  due  to  surface  states  or  collective 
surface  excitations  are  ignored. 

The  expressions  for  the  classical  image  potential  in  an  ideal  metal-insulator-semiconductor  (MIS) 
junction  similar  to  that  shown  in  Fig.  1  diverge  at  the  semiconductor-vacuum  and  vacuum-metal 
interfaces  These  unphysical  divergences  may  be  circumvented  by  employing  an  effective 
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Thomas-Fermi  theory  of  dielectric  screening  to  derive  a  semiclassical  expression  for  the 

image  potential  experienced  by  a  tunneling  electron.  Specifically,  the  Thomas-Fermi  dielectric 
function  in  a  bulk  metal  is  given  by  (£«(^)/ =  d  + /^  / where  is  the  corresponding 
Thomas-Fermi  wave  vector.  Within  the  same  approximation,  the  screening  response  of  an 
intrinsic  bulk  semiconductor  may  be  modeled  by  the  dielectric  function 


{£,{k)l  £^)  = 


1-f- 


\  +  {£j  !  y]) 


(7) 


where  {£^£q)  is  the  long-wavelength  static  dielectric  constant  and  Ys  an  effective  Thomas-Fermi 
wave  vector  which  depends  on  the  total  valence  charge  density.  It  can  then  be  shown  that  the 
resulting  semiclassical  expression  for  the  image  potential  experienced  by  a  point  charge  q  located  at 
position  z  <0  in  the  semiconductor  is 


where 


f: 


1 


(9) 


coincides  with  the  classical  solution  d7.l8] 


represents 
we  define, 


an  additional  response  arising  from  short-range  screening  in  the  semiconductor.  Here, 
in  addition  to  /J  in  (4)  above. 


a„{k) 


k 


(11) 


a.{k)  = 


k 


(12) 


5 


(13) 


l-ajk) 
l  +  a«(^) 


'  l  +  a,(fc) 


When  q  is  located  in  the  vacuum  region  0  <  z  <  s,  the  corresponding  expression  for  the 
semiclassical  image  potential  becomes 


where 


B:  =  -ri„e- 


rif: 


c;  = 


F:  =  {\  +  fi) 


l-j3e‘ 


Finally,  for  q  in  the  metal  (z  >  5 )  one  finds 


l/,(r)  =  exp[(/:^  +  z)]D;rf/:, 


6 


representing  the  effect  of  short-range  screening.  It  is  easily  verified  that  the  formulae  for  Vs  (z), 
V^v  (z),  and  (z)  given  above  reproduce  the  expected  classical  results  in  the  limits  Js^m  and 
that  V'm  (z),  in  panicular,  vanishes. 

An  effective  one-electron  potential  is  obtained  throughout  the  junction  by  adding  the  appropriate 
regional  solution  for  the  semiclassical  image  potential  to  the  band  diagram  shown  in  Fig.  1.  We 
first  consider  the  case  illustrated  in  Fig.  1(b)  when  the  applied  bias  VpB  has  been  chosen  to 
eliminate  any  electric  field  between  tip  and  sample  in  the  absence  of  image  corrections.  Selecting 
the  bottom  of  the  semiconductor  conduction  band  as  an  energy  reference,  the  one-electron  potential 
is  then  given  by  Vs  (z)  in  the  semiconductor,  V^(z)  x  in  the  vacuum  gap,  and  Vm(z)  +  x  •  (^ 
+  Ef)  in  the  metal.  This  function,  though  finite,  is  still  discontinuous  at  the  semiconductor- 
vacuum  and  vacuum-metal  boundaries.  There  are  a  number  of  reasons  for  'his  discontinuity, 
which  originates  from  inadequacies  in  the  semiclassical  treatment  at  distances  within  a  screening 
length  of  the  interfaces.  First,  the  treatment  assumes  bulk  values  for  the  screening  parameters  in 
each  medium.  Second,  the  correct  behavior  of  the  total  effective  potential  (or  self  energy) 
experienced  by  an  electron  moving  from  inside  a  polarizable  medium  to  the  region  immediately 
outside  it,  can  only  be  deduced  by  taking  full  account  of  many-body  effects  Finally,  there  is 
no  a  priori  relationship  between  the  effective  band  bending  produced  by  image  effects  in  the 
immediate  vicinity  of  an  interface,  and  the  one-electron  potential  used  to  model  band  structure, 
which  guarantees  that  the  contribution  from  surface  dipole  layers,  together  with  exchange  and 
correlation  effects,  will  produce  a  smoothly  varying  total  potential  across  the  interfaces.  Since  the 
true  one-electron  potential  is  presumably  continuous  when  all  of  these  effects  are  self-consistently 
accounted  for,  we  supplement  the  semiclassical  picture  by  employing  a  third-order  polynomial  to 
smoothly  join  the  calculated  potential  and  its  derivative  at  points  a  distance  Ysjn^  either  side  of 
the  semiconductor-vacuum  and  vacuum-metal  boundaries.  This  interpolation  is  then  used  to 
represent  the  tunneling  potential  in  the  near-interfacial  regions.  An  example  of  this  interpolated 
potential  is  presented  in  Fig.  2(a)  for  a  junction  consisting  of  a  n-type  hydrogen-terminated  Si(lll) 
surface  and  a  gold  tip  electrode  (with  y;„=1.7A-^  ys=1.16A-l,  and  £5=11.9  1^5-27]^^ 
contrast  our  results  with  the  classical  calculation. 

The  more  general  situation  occurs  when  Vjt  Vfs,  and  there  is  a  non-vanishing  field  in  the  vacuum 
region.  In  the  absence  of  surface  states  this  field  penetrates  the  semiconductor  interior,  and  the 
resultant  band  bending  due  to  the  tip-induced  space  charge  is  schematically  illustrated  in  Fig.  1(a) 
for  n-type  material  when  AO  =  0  -  (x  +  0n)  >  0.  The  surface  potential  qV^i  is  a  sensitive 
function  of  both  gap  spacing  and  bias  voltage,  and  is  given  in  the  depletion  approximation  by 
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(21) 


qV, {s,  V^)  =  A0{[l  +  {sls,)-  +  qVI  Ad))"'  -  ^  /  ^o}'- 

The  resulting  electrostatic  potential  in  the  semiconductor  interior,  v/Tr),  has  finite  range,  vanishing 
at  2  -  -\V,  and  follows  the  well-known  form 

V^(2)  =  V,(l  +  2/^V)^  (22) 

for  -W  <z  <0.  The  depletion  width  W  is  a  function  of  both  the  bulk  donor  density  Nd  and  the 
semiconductor  surface  potential  qV^  through 

W(qV,)  =  {2eyjqNj'\  (23) 

while  the  length  scale  sq  is  obtained  from  the  zero-bias,  zero-separation  depletion  width  according 
to  ^221 


=  (24) 

The  total  one-electron  potential  now  includes  a  contribution  from  y/^z)  as  well  as  that  from  (z) 
for  2  <  0.  As  above,  we  interpolate  within  the  near-interfacial  region  to  obtain  a  continuous  result 
as  a  function  of  position,  and  an  example  calculated  for  the  same  hydrogen-terminated  Si(l  1 1)/Au 
junction,  at  a  sample  bias  of  -100  mV  with  =  5  x  10  (2  cm'^,  is  illustrated  in  Fig.  2(b),  where 
the  additional  long-range  effect  of  Debye  screening  due  to  free  carriers  has  been  ignored. 

To  compute  a  tunneling  current  density,  the  appropriate  potential  is  substituted  into  a  one¬ 
dimensional  effective-mass  Schrodinger  equation  for  electrons  originating  from  the  conduction 
band  of  the  semiconductor.  This  equation  is  then  solved  numerically  to  evaluate  a  transmission 
probability  D{Ex,  s),  where  Ez  is  the  energy  component  normal  to  the  barrier.  The  net  tunneling 
current  density  is  then  given  by  the  expression 

J(s)  =  -^fdE,D{E,,s)[N,(E,)-N,{E,-qV)]  .  (25) 

y  ■'0 

where  Nq  is  the  one-dimensional  supply  function  for  electrons  at  finite  temperature  in  the 
semiconductor  and  metal,  respectively,  while  y  represents  an  anisotropy  factor  for  the  constant 
energy  surfaces  of  the  six  equivalent  Si<l(X)>  ellipsoidal  pockets  at  the  indirect  conduction  band 
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minimum  projected  along  the  (111)  direction.  The  resulting  effective  mass  for  tunneling  in  the 
(111)  direction  is  mz*  =  '{mi  =  0.26mo.  Further  details  concerning  the  calculation  of  DiEz,  s)  may 
be  found  in  Ref.  [22]. 

III.  Results 

The  tunneling  current  density  as  a  function  of  gap  spacing  calculated  on  the  basis  of  Eq.  (25), 
using  the  one-electron  effective  potential  of  Fig.  2(b),  is  presented  in  Fig.  3  for  a  representative 
bias  voltage  which  produces  substantial  majority  carrier  depletion  in  the  near  surface  region  of  the 
semiconductor.  For  comparison,  the  results  obtained  using  the  potential  of  Fig.  1(a),  neglecting 
image  effects,  are  also  shown.  At  each  of  the  doping  levels  considered,  we  see  roughly  an  order 
of  magnitude  increase  in  the  overall  current  density  at  any  reasonable  gap  separation  due  to  the 
inclusion  of  image  forces.  The  shape  of  the  j  -  s  curves  on  a  semi- logarithmic  plot  deviates  at 
small  s  from  the  straight  line  predicted  for  vacuum  tunneling  because  of  the  influence  of  tip- 
induced  space  charge  on  current  transport  through  the  junction,  as  explained  previously  What 
we  wish  to  call  attention  to  here  is  that  the  j  -  s  characteristics  for  a  given  doping  display  a  similar 
shape  whether  or  not  image  effects  are  included. 

This  observation  is  more  quantitatively  expressed  in  Fig.  4,  where  an  apparent  barrier  height  as  a 
function  of  gap  spacing  is  obtained  via  Eq.  (1)  from  the  curves  in  Fig.  3.  We  see  that  the  barrier 
height  displays  nearly  the  same  behavior  as  a  function  of  s  (for  a  given  doping)  with  or  without 
image  forces,  but  that  including  image  effects  results  in  an  overall  lowering  of  the  apparent  barrier 
by  a  few  tenths  of  an  eV.  Also  presented  in  Fig.  4  are  the  doping-independent  results  for  0a 
when  the  applied  bias  has  been  chosen  to  eliminate  any  tip-induced  band  bending  as  shown  in 
Fig.  1(b).  It  is  clear  from  the  figure  that  the  principal  influence  on  0a  arises  from  tip-induced 
space  charge,  and  that  the  image  potential  constitutes  a  relatively  small  perturbation. 

The  influence  of  image  forces  on  the  tunneling  current- voltage  characteristics  at  constant  separation 
may  be  calculated  in  a  manner  similar  to  the  j  -  s  characteristics  presented  above,  and  the  results  are 
shown  in  Fig.  5,  as  a  function  of  doping,  for  s  =  5A.  To  facilitate  a  comparison  between  j  -  V 
curves  calculated  with  and  without  image  corrections  at  fixed  doping,  and  between  j  -V  curves 
obtained  for  different  doping  levels,  both  the  current  density  and  applied  bias  have  been  scaled  to 
their  respective  flat  band  values.  As  already  noted,  the  image  potential  increases  the  current  density 
for  a  given  applied  bias  and  electrode  separation,  but  we  see  here  that  the  overall  shape  of  the  j  -V 
characteristic  remains  essentially  unchanged.  This  shape  is,  again,  dominated  by  the  tip-induced 
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space  charge.  A  discussion  of  the  somewhat  unusual  properties  of  the  reverse-bias  characteristics 
may  be  found  in  the  literature 

IV.  Discussion 

Our  results  for  the  specific  case  of  a  planar,  passivated  n-type  Si  (1 1 1)/Au  junction  indicate  that  the 
influence  of  image  effects  on  the  apparent  barrier  height  is  small,  and  may  not  be  experimentally 
detectable  in  practice.  The  role  of  tip-induced  space  charge  in  determining  the  apparent  barrier 
decreases  with  increasing  bias  as  flat  band  conditions  are  approached  so  that  image  effects,  if 
detectable,  will  be  most  prominent  at  bias  voltages  very  near  The  shape  of  the  tunneling 
current-voltage  characteristics  also  appears  to  be  insensitive  to  image  forces. 

In  what  follows,  we  examine  these  conclusions  from  the  perspective  of  a  WKB  analysis  of 
tunneling  through  the  compound  barrier  created  either  by  the  space  charge  in  the  semiconductor 
neglecting  image  effects,  followed  by  the  vacuum  gap,  or  the  image  potential  in  the  semiconductor 
neglecting  space  charge,  and  the  vacuum  gap.  In  this  way  we  can  identify  important  qualitative 
features  of  the  realistic  junction  problem  containing  both  elements  which  will  give  us  some  insight 
into  the  relevance  of  specific  material  parameters  to  our  results,  and  the  circumstances  under  which 
image  effects  might  play  a  more  significant  role. 

The  simplest  case  arises  when  the  only  barrier  which  need  be  considered  is  due  to  the  vacuum  gap. 
If  we  denote  by  D''3C(£^,  v)  the  energy-dependent  transmission  coefficient,  then  in  the  WKB 
approximation  one  obtains 


lnD''*'(E,,5)  =  -2(2/n,  /  n^f\x  " 


(26) 


The  differential  current  density  per  unit  energy  at  energy  Ez  is  then 

dj{E,,s)ldE,=D''*\E,,s)[N,{E,)-N,{E,-qV)\,  (27) 


and  the  condition  for  an  extremum  in  this  quantity  is  simply 

-m  kT)  +  d\nD'''‘{E„s)l  dE\  =0- 


(28) 
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This  condition,  when  appropriated  generalized,  gives  Eq.  (29)  and  is  relevant  to  the  discussion 
following  (31).  However,  in  the  present  case,  the  extremum  in  (28)  turns  out  to  be  a  local 
minimum,  and  a  maximum  is  found  at  the  end  point  £„  =  0.  If  the  differential  current  density  at 
this  maximum  is  used  to  characterize  the  integrated  mnneling  current,  Eq.  (2)  yields  fpA  =  X  for  the 
apparent  barrier  height. 

The  next  situation  which  is  straightforward  to  analyze  arises  when  we  consider  the  effect  of 
semiconductor  space  charge  and  ignore  the  vacuum  barrier.  This  corresponds  to  the  problem  of 
tunneling  through  an  ideal  Schottky  barrier  structure.  In  this  case,  a  WKB  analysis  reveals  that  the 
peak  of  the  differential  current  distribution  occurs  at  an  energy 

EJV,={cosh[q{h^NJAm,e,f^/kT]Y\  (29) 

In  the  limit  of  low  doping  and  long  depletion  widths,  this  ratio  approaches  unity,  so  that  the 

IkT 

differential  current  density  is  proponional  to  ^  '  and  one  recovers  the  thermionic  emission 
picture  of  current  transport  across  a  Schottky  barrier 

The  combined  effect  of  the  semiconductor  space  charge  and  a  vacuum  barrier  is  more  difficult  to 
follow  through  completely,  but  the  qualitative  features  are  apparent  in  at  least  one  important  limit. 
As  noted  in  connection  with  Eq.  (21),  the  semiconductor  surface  potential  is  a  function  of  the 
vacuum  gap  width  s  .  If  one  considers  the  physically  realistic  circumstance  of  sIsq  «  1,  which 
corresponds  to  electrode  separations  that  are  small  compared  to  the  sample  depletion  width,  the 
energy  distribution  will  be  almost  completely  determined  by  the  semiconductor  space  charge.  A 
specific  example,  based  on  the  more  accurate  calculations  described  in  Section  II,  is  presented  in 
Fig.  6,  where  we  see  that  Em/Vd  is  nearly  one,  independent  of  s.  Thus,  for  low  doping,  current 
transport  is  once  again  well-described  by  a  thermionic  emission  picture,  and  one  finds  from  Eq, 
(2)  that 


=  ^1/2  ^  !  kT)[dV,  1 ,95].  (30) 

Since  the  surface  potential  decreases  with  increasing  separation,  r,  the  second  term  in  (30)  is 
opposite  in  sign  to  the  vacuum  barrier  contribution.  The  distance  dependence  of  Vd  then  has  a 
profound  influence  on  the  apparent  barrier  height  through  its  effect  on  the  tunneling  energy 
distribution. 
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Finally,  we  turn  to  the  question  of  image  forces  and  consider  the  case  where  we  can  neglect  tip- 
induced  band  bending.  A  key  issue  is  the  extent  to  which  image  effects  in  the  semiconductor 
interior  perturb  the  energy  distribution  of  tunneling  electrons  relative  to  that  characteristic  of  the 
vacuum  barrier  alone.  Ignoring  image  corrections  in  the  vacuum  region  for  simplicity,  the 
extremum  condition  of  Eq.  (28)  now  becomes. 


-{\l  kT)  +  d\nD''^"{E„s)l  dE^  +^lnD‘™^£„5)/^£, 


=  0 


(31) 


which  includes  a  contribution  from  the  image-induced  barrier  in  the  semiconductor.  A  useful 
estimate  of  the  last  term  in  (3 1)  is  provided  by  considering  the  limit  of  infinite  separation  between 
tip  and  sample.  The  semiconductor  barrier  is  then  a  screened  coulomb  potential  due  to  a  single 
image  charge,  as  appropriate  for  the  free  semiconductor- vacuum  interface.  Since  the  proximity  of 
the  metal  counter-electrode  creates  multiple  image  charges  whose  net  effect  limits  the  range  of  the 
semiconductor  barrier,  we  may  reasonably  expect  this  approximation  to  provide  an  upper  bound  to 
the  image-induced  perturbation  of  the  tunneling  energy  distribution.  Employing  the  well-known 
expression  for  the  energy-dependent  WKB  transmission  coefficient  through  a  coulomb  barrier 
the  relative  significance  of  semiconductor  and  vacuum  terms  may  be  assessed  via  the 
dimensionless  ratio 


(  •  ^ 

m 

r  ^ 

jUJ 

(32) 


For  typical  values  of  the  electron  affinity  in  semiconductors,  the  last  term  in  parenthesis  is  of  order 
(1  A/s),  while  the  maximum  possible  value  for  the  remaining  parameters  is  0.26.  Evaluation  for  the 
specific  case  of  Si(l  1 1)  gives  0.05  (lA/s)  for  this  ratio.  Thus,  at  gap  separations  of  order  ten 
angstroms  or  more,  the  corrections  to  the  maximum  of  the  energy  distribution  are  expected  to  be 
small  for  any  material,  and  certainly  negligible  for  silicon.  We  can  inquire  further,  under  the 
presumption  that  the  semiconductor  barrier  will  not  be  negligible  under  all  circumstances,  what 
may  at  most  be  expected  for  the  shift  in  the  location  of  the  energy  maximum.  As  an  extreme 
circumstance  we  retain  the  screened  coulomb  potential  in  the  semiconductor,  while  completely 
ignoring  the  vacuum  barrier,  and  find 
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For  Si  (1 1 1)  at  room  temperature,  [EmlkT)  =  0.76,  whereas  for  optimum  parameters  it  will  be  no 
more  than  a  factor  of  3  greater. 

The  above  considerations  suggest  therefore,  that,  relative  to  ^  zero-energy  WKB  analysis  is  an 
appropriate  starting  point  for  discussing  how  image  forces  influence  the  apparent  barrier  height  in 
this  system.  That  supposition  is  supported  by  numerical  results  obtained  from  the  detailed 
formalism  presented  in  Section  II.  In  particular,  Fig.  7  displays  the  differential  current  density  per 
unit  energy,  as  a  function  of  gap  spacing,  calculated  for  the  semiclassical  multiple-image  interaction 
using  the  one  electron  potential  of  Fig.  2(a).  As  one  would  expect,  these  distributions  show 
resonant  transmission  for  electron  energies  very  near  the  conduction  band  minimum,  due  to  the 
shallow,  long-range  attractive  ponion  of  the  image  potential  in  the  semiconductor  interior,  and, 
funhermore,  they  peak  just  below  (EzikT)  =  1,  as  inferred  above.  Of  particular  note  in  Fig.  7  is 
that,  in  marked  contrast  to  the  situation  for  tip  induced  band  bending  (Fig.  6),  the  energy 
distribution  for  tunneling  electrons  appears  to  be  only  weakly  dependent  on  s  (at  least  when  the 
semiconductor  barrier  is  of  finite  range),  a  relevant  point  outside  the  scope  of  our  arguments  thus 
far. 

As  a  consequence  of  the  scaling  property  (5)  for  the  classical  image  potential  in  a  semiconductor, 
the  zero-energy  WKB  transmission  coefficient  through  this  barrier  may  be  parameterized  in  terms 
of  the  dimensionless  integral 


(34) 


It  follows  that  the  contribution  of  the  semiconductor  to  the  apparent  barrier  height  at  =  0  is  given 
by 


<9lnD”""°"\£,  =0,5)  _ 


ds 


=  f(/J) 


m 


lm„ 


^Ke,s  j 


1/2 


(35) 


In  relation  to  the  vacuum  barrier  then,  this  represents  a  correction  of  order 


d\nD''*\E,=^Q,s)l  ds 


\^0  J 


,2  \ 
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(36) 
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Using  our  estimate  above,  the  last  term  in  parentheses  is  roughly  Now  1(13)  has  a 

maximum  value  of  0.51  at  /3  =1  (an  MVM  junction),  and  decreases  smoothly  toward  zero  in  the 
lim.it  that  P  vanishes.  For  any  reasonable  choice  of  the  semiconductor  dielectric  constant, 
however,  /( p)  e.xceeds  0.4.  Thus,  in  general,  our  correction  term  can  at  most  be  of  order 
0. 3(1  A/s) and  for  Si  (1 1 1)  is  only  a  factor  of  two  smaller.  Over  a  physically  reasonable  range 
of  gap  spacings  from  10  -  1  A,  for  example,  the  effect  of  the  semiconductor  image  potential,  in  this 
approximation,  will  be  to  increase  the  apparent  barrier  height  between  5  and  15%. 

Because  of  our  reliance  on  a  zero-energy  analysis,  the  above  estimate  presents  an  upper  limit  to  the 
magnitude  of  the  effect  one  may  reasonably  expect.  Indeed,  our  numerical  calculations  for  Si(l  1 1) 
show  that  the  inclusion  of  image  forces  produces  a  small  net  decrease  in  the  apparent  barrier  height 
due  to  a  reduction  of  the  vacuum  barrier,  which  is  fully  consistent  with  the  results  for  MVM 
tunneling. 


One  may  pause  to  consider  at  this  point  whether  the  situation  will  be  fundamentally  altered  if  we 
replace  the  far  metal  electrode  with  a  second  semiconductor  whose  polarizability  is  less  than  or 
equal  to  that  of  the  first.  In  this  case,  the  integral  (34)  diverges  (with  *”)’  since  there  is  no 

longer  any  change  in  sign  for  U^p,  P',q)  and  the  image-induced  barrier  is  long  range.  Reflecting 
on  (33),  however,  we  are  reminded  that  this  barrier  will  be  probed  not  at  =  0,  but  rather  at  = 
kT.  The  effective  range  of  the  image  interaction  at  this  energy  is  then  of  order 


\6Ke^kT  J 


(37) 


which,  at  r  =  300  K,  is  less  than  25  A  for  optimum  p.  Thus,  there  is  nothing  in  our  arguments 
which  leads  us  to  suspect  circumstances  will  be  dramatically  different  for  semiconductor-vacuum- 
semiconductor  (SVS)  tunneling,  unless  one  can  arrange  to  conduct  the  experiments  at  low 
temperature. 

V.  Summary  and  Conclusions 

We  have  investigated  the  influence  of  the  image  potential  on  tunneling  through  a  planar  metal- 
vacuum-semiconductor  Junction.  This  problem  differs  qualitatively  from  the  metal-vacuum-metal 
junction  because  there  is  now  an  image  potential  in  the  semiconductor  interior  as  well  as  in  the 
vacuum  region.  The  classical  multiple-image  expression  for  an  MVS  system  has  been  extended. 


using  a  Thomas-Fermi  like  dielectric  function  in  the  semiconductor  together  with  the  usual 
Thomas-Fermi  theory  for  the  metal,  to  provide  expressions  which  are  everywhere  finite.  The 
resulting  semiclassical  image  potentials  in  semiconductor,  vacuum,  and  metal  must  still  be 
smoothly  joined,  however,  by  an  interpolation  through  the  interfaces.  Using  these  potentials,  the 
tunneling  current  at  finite  temperature  was  calculated  as  a  function  of  separation  for  the  panicular 
case  of  a  n-type  hydrogen-passr  •’ted  Si(l  11)  -  Au  system,  and  an  apparent  barrier  height  0^  then 
extracted  to  simulate  the  results  of  STM  modulation  experiments. 

To  obtain  a  more  general  understanding  of  these  detailed  calculations,  qualitative  arguments  based 
on  a  WKB  picture  have  also  been  presented.  If  the  image  potential  is  neglected,  <X>^  is  well 
described  by  terms  involving,  first,  the  semiconductor's  electron  affinity,  and,  second,  the  effect 
of  tip-induced  band  bending.  When  the  contribution  of  the  image  potential  within  the 
semiconductor  is  separately  considered,  a  WKB  estimate  indicates  it  will  perturb  the  apparent 
barrier  height  by  at  most  10%  at  reasonable  tip-sample  distances.  This  crude  estimate  helps  one 
understand  the  principal  result  of  the  detailed  calculations  presented  in  Figures  3  and  4:  Although 
the  tunneling  barrier  for  an  MVS  junction  differs  from  that  for  an  MVM  junction,  the  image 
potential  in  the  semiconductor  makes  only  a  modest  contribution  to  the  apparent  barrier  height. 
Funhermore,  as  our  calculations  reveal,  the  dominant  effect  is  a  slight  lowering  of  0a  .  due  to  the 
attractive  image  potential  in  the  vacuum  region  in  conformity  w  ith  expectations  based  on  MVM 
tunneling,  rather  than  any  increase  in  0a  from  the  repulsive  image  potential  within  uie 
semiconductor. 

One  may  have  thought  that  the  semiconductor's  internal  image  potential,  being  long-range,  would 
have  a  substantial  effect  on  carriers  originating  from  the  bulk  conduction  band  edge,  and  thus 
influence  not  only  tunneling  currents,  but  the  apparent  barrier  height  as  well.  Our  calculations, 
together  with  the  WKB  analysis,  indicate  otherwise.  Three  qualitative  arguments  account  for  this 
observation;  First,  in  the  case  of  flat  bands,  the  potential  due  to  polarization  charge  at  the  vacuum- 
metal  interface  competes  with  that  of  the  semiconductor-vacuum  interface  to  produce  a  barrier  in 
the  semiconductor  interior  which  is  short  range;  Second,  for  the  more  general  case  which  includes 
a  semiconductor  diffusion  potential  in  either  depletion  or  weak  inversion,  the  overwhelming 
influence  on  the  carrier  distribution  arises  from  the  diffusion  potential,  and  not  the  image  potential. 
In  the  case  this  diffusion  potential  is  tip-induced,  we  have  shown  explicitly  that  0a  is  almost 
completely  dominated  by  the  distance  dependence  of  field  penetration  in  the  semiconductor,  and 
that  image  effects  are  of  relatively  minor  importance;  Finally,  for  flat  band  conditions,  most  of  the 
tunneling  current  is  provided  by  electrons  whose  energy  is  of  order  kT  above  the  band  edge.  For 
these  electrons,  the  image-induced  barrier  is  effectively  shon  range  at  room  temperature,  even  if 
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one  neglects  the  influence  of  polarization  charge  at  the  metal  tip.  The  same  argument  applies,  in 
principle,  to  the  semiconductor-vacuum-semiconductor  junction,  but  whether  or  not  the  image 
potential  will  have  a  significant  influence  in  this  system  at  low  temperatures  remains  an  open 
question. 
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Figure  2: 


Figure  3: 


Figure  4: 


Figure  5: 


Figure  6: 


Figure  7: 


Figure  Captions 


(a)  Schematic  energy  diagram  for  an  ideal  semiconductor-vacuum-metal  junction  at 
zero  applied  bias  V,  neglecting  image  effects.  0  is  the  metal  work  function,  Ep  the 
metal  Fermi  energy,  x  the  semiconductor  electron  affinity,  and  (pn  the  offset  from  the 
Fermi  level  to  the  tottom  of  the  semiconductor  conduction  band.  Tlie  drawing  is  not  to 
scale  since  the  depletion  width  W  is  typically  much  greater  than  the  gap  spacing  s,  and 
X  exceeds  the  band  gap,  Ec-  Ev . 

(b)  Flat  band  situation  when  the  applied  bias  VfB  has  been  chosen  to  eliminate  th*’ 

electric  field  arising  from  a  difference  in  metal  and  semiconductor  work  functions, 

<P-(  <Pn)- 


(a)  Effective  one-electron  potential  energy  as  a  function  of  position  for  V  =  VfB  and 

s  =  10  A.  Solid  line:  including  the  semiclassical  image  potential  with  interpolation. 
Dotted  line:  including  the  classical  image  interaction.  Dashed  line:  neglecting  image 
effects.  X  -  ^  The  potentials  are  independent 

of  doping. 

(b)  Effective  one-electron  potential  energy  as  a  function  of  position  at  a  representative 
bias  which  produces  depletion  in  the  semiconductor  =  5  x  10^^  cm'^).  All  other 
parameters  as  above.  Solid  line:  including  the  semiclassical  image  potential  with 
interpolation.  Dashed  line:  neglecting  image  effects.  Note  the  offset  of  the 
semiconductor  surface  potential  from  the  conduction  band  edge  in  the  bulk. 


Tunneling  current  density  j  versus  electrode  separation  s  as  a  function  of  doping  for 
fixed  bias  and  T  =  300  K.  Solid  line:  including  the  semiclassical  image  potential  as  in 
Figure  2(b).  Dashed  line:  neglecting  image  effects. 


Apparent  barrier  height  0^  versus  electrode  separation  r  as  a  function  of  doping.  Solid 
line:  including  the  semiclassical  image  potential.  Dashed  line:  neglecting  image 
effects.  The  flat  band  results  are  essentially  independent  of  doping  in  the  non¬ 
degenerate  regime. 


Tunneling  current  density  j  versus  applied  bias  as  a  function  of  aoping  for  fixed 
separation  and  T  =  300  K.  Solid  line:  including  the  semiclassical  image  potential. 
Dashed  line:  neglecting  image  effects.  The  current  density  and  applied  bias  are  scaled 
to  their  respective  flat  band  values  to  facilitate  a  comparison  of  the  shape  of  the  ;  - 
characteristics. 


Differential  tunneling  current  density  distribution  versus  normal  energy  Ez  (scaled  to  the 
semiconductor  surface  potential)  as  a  function  of  electrode  separation,  neglecting  image 
effects,  for  V  =  -100  mV  and  =  5  x  10^^  cm'^.  Inset:  tip-induced  semiconductor 
surface  potential  versus  separation  at  the  same  fixed  bias. 


Differential  tunneling  current  density  distribution  versus  normal  energy  Ez  (scaled  to 
kT)  as  a  function  of  electrode  separation,  including  the  semiclassical  image  potential  in 
the  absence  of  tip-induced  band  bending  iV=  VfB). 
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Section  II 

INSTRUMENTATION 


Figure  Captions 


Figure  l(a*b):  UHV  STM  system  in  Surface  Physics  Laboratory  at  Texas  A&M 

University.  The  entire  assembly  rests  atop  a  specially  engineered  vibration  isolation  platform. 

Figure  2(a):  Orthographic  views  of  main  analytical  chamber.  Ranges  are  arranged  along 

the  chamber  axis  in  three  clusters.  From  left  to  right:  Surface  analytical  cluster  with  provision  for 
Auger  Spectroscopy,  LEED,  X-ray  Photoelectron  Spectroscopy,  and  Ion  Sputtering;  Sample  pick- 
off  from  preparation  chamber  /  entry  lock  with  provision  for  sample  storage;  STM  cluster. 
Maximum  flexibility  is  maintained  by  the  orthogonal  arrangement  of  8"  ports  in  the  STM  cluster, 
permitting  bottom,  top,  or  side  mounting  of  an  STM  flange.  Introduction  of  samples  to  the  main 
chamber  and  movement  along  the  axis  of  this  chamber  is  via  orthogonal,  magnetically  coupled 
linear/rotary  transporters. 

Figure  2(b):  Schedule  describing  location,  orientation,  and  function  of  various  ports 

available  on  the  analytical  chamber,  as  well  as  sample  target  positions. 

Figure  3(a):  Orthographic  views  of  UHV  pumping  assembly  consisting  of  Perkin-Elmer 

220  LAsec  D-I  ion  pump,  titanium  sublimation  pump  with  liquid  nitrogen  cooled  cryoshroud,  and 
custom  mounting  for  SAES  nonevaporable  getter  (NEG)  module.  Auxiliary  flanges  for  ion  gauge 
and/or  quadrupole  mass  spectrometer. 

Figure  3(b):  Small  volume  sample  preparation  chamber  for  directed  gas-phase  dosing, 

electrochemical  deposition  cells,  or  filament  evaporators.  A  gate  valve  separating  this  side  chamber 
from  the  main  analytical  volume  permits  dosing  with  high  vapor  pressure  substances  without  the 
ri.sk  of  either  contaminating  analytical  instruments  or  exposing  the  STM  tip.  A  30  L/sec  Perkin- 
Elmer  D-I  ion  pump  and  Balzer's  60  L/sec  turbo-molecular  pump,  attached  to  the  entry  lock, 
provide  pumping  for  this  isolated  segment. 

Figure  4(a):  Ambient  (vertical)  vibration  amplitude  versus  frequency,  on  laboratory  floor 

(upper  trace)  and  atop  vibration  isolation  platform  (lower  trace).  Vertical  scale:  logarithmic,  from  - 
60  db  (=  1  pm)  to  -t-4()  db  (=  100  nm),  with  10  db  per  division.  Dashed  line  indicates  0  db  (=  1 
nm)  displacement.  Horizontal  scale:  linear,  from  0  Hz  to  50  Hz,  with  5  Hz  per  division  and  a 
frequency  resolution  of  0.125  Hz.  Platform  resonance  occurs  at  1.25  Hz  and  unity  gain  crossover 
in  the  transfer  function  occurs  at  1.625  Hz.  The  1  nm  displacement  crossover  frequency  is  4.00 
Hz.  Above  1 .625  Hz,  the  platform  successfully  filters  out  vertical  seismic  disturbances,  reaching 
an  ultimate  suppression  factor  greater  than  40  db  beyond  10  Hz. 

Figure  4(b):  Ambient  (horizontal)  vibration  amplitude  versus  frequency,  on  laboratory 

floor  (upper  trace)  and  atop  vibration  isolation  platform  (lower  trace).  Vertical  and  horizontal 
scales  as  in  4(a).  Platform  resonance  occurs  at  1. 125  Hz  and  unity  gain  crossover  in  the  transfer 
function  occurs  at  1.50  Hz.  The  1  nm  displacement  crossover  frequency  is  3.00  Hz  in  this  case. 
Ultimate  suppression  above  10  Hz  is  between  30  db  and  40  db,  depending  on  frequency. 

Figure  5(a):  Finite-element  simulation  (ALGOR)  of  normal  modes  for  magnetically 

coupled  linear/rotary  transporters  fixed  at  one  end  (cantilevered  support).  Lowest  vibrational 
resonance  occurs  at  approximately  25  Hz. 

Figure  5(b):  Finite-element  simulation  of  normal  modes  for  magnetically  coupled 

linear/rotary  transporters  fixed  at  both  ends.  Lowest  vibrational  resonance  has  been  shifted  to 
approximately  87  Hz. 


Figure  Captions  (cont'd) 


Figure  6(a-b):  Normal  modes  of  center-  and  transverse-axis  magnetically  coupled 

linear/rotary  transporters  in  rigidized  support  scheme  implemented  in  Figure  1.  Note  that  lowest 
resonant  frequency  is  close  to  that  of  ideal  fixed-fixed  boundary  conditions  in  Figure  5(b)  above. 

Figure  7(a-b):  Close-up  view  of  Omicron  STM  which  shows  flange  mounting,  single- 

stage  spring  suspension  with  eddy-current  damping,  coarse  two-dimensional  sample  positioner 
(left  center),  piezo-tripod  with  exchangeable  tip  holder  (right  center),  and  BEEM  current 
preamplifier  (bottom  center).  Hidden  from  view  by  supports  for  spring  suspension  is  tunneling 
current  preamplifier  mounted  directly  behind  the  piezo  scanner. 

Figure  8:  View  along  main  axis  of  analytical  chamber,  through  port  1 ;  Looking  past 

Cylindrical  Mirror  Analyzer  (in  foreground),  infra-red  bake  out  lamps,  and  sample  carousel, 
toward  STM  (in  background)  with  stage  at  transfer  elevation,  and  pincer-grip  wobble  stick. 
Center-axis  transponer,  holding  tip  or  sample,  is  positioned  ju.st  above  STM. 

Figure  9(a-b):  Cross-section  through  main  chamber  at  pons  6,7,8  and  9,  illustrating 

proposed  pick-off  scheme  after  sample  introduction  from  entry  lock  or  preparation  chamber. 
Wobble  stick,  sample  plattens,  and  carousel  are  Omicron  design.  Wobble  stick  includes  pincer 
movement  for  positive  grip  on  platten.  Sample  plattens  are  shown  mounted  on  custom  rack 
attached  to  the  end  of  a  magnetic  linear/rotary  transporter  with  30”  stroke.  Overhead  view  port 
allows  convenient  inspection  of  transfer  process  and  loaded  carousel.  Wobble  stick  retraction  is 
sufficient  for  plattens  to  be  exchanged  from  front  to  back  on  carousel,  thereby  permitting  a  reversal 
of  their  orientation  relative  to  chamber  center  line  once  the  carousel  is  rotated.  This  is  crucial  as 
.sample  and  tip  loading  must  proceed  with  opposite  orientation  in  the  STM.  Rotation  of  carousel  is 
accomplished  by  tapping  a  clutch  lever,  which  is  inside  a  20  degree  cone  angle,  with  wobble  stick. 
A  high  conductance,  stainless  steel  me.sh  screen  serves  to  protect  the  pumping  assembly  from  any 
inadvertently  dropped  plattens. 

Figure  10(a):  Cross-section  through  main  chamber  at  ports  10,11,12  and  13,  showing 

relative  orientation  of  center-line  magnetic  linear/rotary  transponer  and  a  second  pincer-grip  wobble 
stick  during  the  initial  stage  of  sample  or  tip  transfer  to  STM.  Overhead  and  side  view  pons 
facilitate  inspection  of  transfer  process. 

Figure  10(b):  Once  magnetic  transporter  has  been  fully  withdrawn,  sample  or  tip  may  be 

loaded  onto  the  STM  stage  by  levelling  wobble  stick  as  shown. 
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Figure  10(a) 
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Section  III 


AUGER  SPECTROSCOPY 


Figure  Captions 


Figure  1:  Auger  spectrum  from  copper  test  sample.  The  sample  was  sequentially  cleaned 

using  trichloroethane,  acetone,  and  then  ethanol,  in  an  ultrasonic  cleaner.  The  spectrum  was  taken 
with  a  Perkin  Elmer  15-155  Cylindrical  Mirror  Analyzer. 


Figure  2:  Auger  spectrum  from  copper  test  sample.  This  part  of  the  sample  was  etched  in 

FeCU  to  remove  oxygen.  Compare  the  ratio  of  peaks  4  and  8  to  the  same  ratio  in  figure  1.  Also 
note  me  ratio  of  peaks  3  and  8  in  each  figure.  The  magnitude  of  peak  2  in  Figure  2  is  due  to  the 
residue  from  the  FeCl^. 
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Section  IV 

TUNNELING  MICROSCOPY 


Figure  Captions 


Figure  1(a):  UHV  test  data  of  a  Si(lll)  7x7  surface  obtained  with  the  instrument 

presently  at  Texas  A&M.  The  sample  was  biased  at  2.0V  with  a  feedback  current  of  0.06nA.  An 
atomic  step  is  visible  at  the  top  of  the  image.  This  data  was  acquired  at  the  Omicron  factory. 

Figure  1(b):  A  grey  scale  rendition  of  the  same  silicon  image. 

Figure  2(a):  Topographic  image  of  highly  oriented  pyrolytic  graphite  (HOPG)  showing 

atomic  resolution  over  a  long  (1(X)  Angstrom)  scan.  The  image  was  obtained  in  air  with  a  platinum 
tip  cut  from  wire.  The  sample  bias  was  50m V  with  a  feedback  current  of  l.OnA. 

Figure  2(b):  Current  image  taken  concurrently  with  the  topographic  image  of  Figure 

2(a).  The  total  time  to  take  the  data  was  1 14  sec;  the  feedback  time  constant  was  set  such  that  the 
atomic  corrugation  modulated  the  current. 

Figure  3(a):  Large  area  scan  (2()0()  Angstroms  by  2000  Angstroms)  of  the  HOPG 

surface  showing  mono-atomic  steps.  The  measured  step  height  is  about  3.4  Angstroms,  in  good 
agreement  with  the  accepted  value  of  3.35  Angstroms. 

Figure  3(b):  A  grey  scale  rendition  of  the  same  HOPG  image. 

Figure  4:  Grey  scale  image  of  a  tip  crash  site  on  HOPG.  This  distinctive  feature 

could  be  imaged  for  hours  and  provided  an  opportunity  to  practice  maneuvering  the  STM  image 
frame. 

Figure  5(a):  A  typical  grey  scale  image  of  the  surface  of  bulk  gold.  The  large,  rounded 

features  such  as  these  are  commonly  observed.  This  image  was  taken  with  a  sample  bias  of 
lOOmV  and  a  feedback  current  of  1  .OnA. 

Figure  S(b):  Another  grey  scale  image  of  bulk  gold.  The  data  was  obtained  in  air  with  a 

sample  bias  voltage  of  2.5V  and  a  feedback  current  of  ().()5nA,  yielding  a  nominal  tunneling 
resistance  of  50  gigaohms.  Large,  rounded  features  are  resolved,  but  finer  details  are  obscured 
presumably  due  to  larger  tip-sample  separation  than  in  5(a). 

Figure  6(a-c):  Typical  current-voltage  (I-V)  characteristics  (left  side  of  each  figure)  for 

junctions  formed  between  a  platinum  STM  tip  and  various  samples:  (a)  graphite,  (b)  bulk  gold, 
and(c)  oxidized  tungsten.  Note  that  (a)  and  (b)  have  qualitatively  different  behavior  near  zero  bias. 
The  characteristics  for  graphite  are  less  linear  than  the  characteristics  for  gold.  This  is  consistent 
with  the  greater  difference  in  the  work  function  values  between  graphite  and  platinum  compared  to 
gold  and  platinum.  The  oxidized  tungsten  characteristics  behave  similarly  to  those  of  a 
.semiconductor.  The  right  side  of  the  figures  shows  the  junction  conductance,  calculated  by 
numerically  differentiating  the  measured  I-V  characteristics.  The  effect  of  noise  in  this  coarse  data 
is  readily  apparent. 


Figure  1(a) 


Figure  1(b) 


Figure  2(a) 


